Luttinger liquid superlattices: realization of gapless insulating phases 
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We investigate LL superlattices, a periodic structure composed of two kinds of one-dimensional 
systems of interacting electrons. We calculate several properties of the low-energy sector: the ef- 
fective charge and spin velocities, the compressibility, various correlation functions, the Landauer 
conductance and the Drude weight. The low-energy properties are subsumed into effective param- 
eters, much like homogeneous one-dimensional systems. A generic result is the weighted average 
nature of these parameters, in proportion to the spatial extent of the underlying subunits, pointing 
to the possibility of "engineered" structures. As a specific realization, we consider a one-dimensional 
Hubbard superlattice, which consists of a periodic arrangement of two long Hubbard chains with 
different coupling constants and different hopping amplitudes. This system exhibits a rich phase 
diagram with several phases, both metallic and insulating. We have found that gapless insulating 
phases are present over a wide range of parameters. 

PACS numbers: 71.10.Pm, 71.10.Fd, 71..30.+h, 73.22.-f, 73.6.3.-b 



I. INTRODUCTION 

The physics of one-dimensional electronic systems has 
been the subject of a vigorous onslaught recently, both 
theoretical and experimental. Experimentally, thft-alDility 
to grow nanostrup|UMjeft-|»iich as quantum wirestlaa and 
carbon nanotubesoQ'Q'QEla has enabled, for the first time, 
the investigation of systems of a truly one-dimensional 
nature. On the theoretical side, the peculiarities of the 
behavior of interacting electrons in one dimension have 
culminated in the proposal p£|a-.miimm^uH^ class 
dubbed the Luttinger HquidBe-BByeO (LL), which 
stands in sharp contrast with the higher dimensional 
Fermi liquids established by Landau. The LL is charac- 
terized by the absence of stable quasi-particles, its low- 
energy sector being exhausted by collective charge and 
spin density excitations. Since the latter travel at dif- 
ferent velocities, an added electron splits up into well 
separated charge and spin degrees of freedom. Further- 
more, correlation functions decay in a power law fash- 
ion, with exponents set by only a few parameters. This 
generic behavior has been tested and confirmed in the 
case of edge transport in sA^stemSpwhirh exhibit the frac- 
tional quantum Hall effect .OLataHJGil LL theory has also 
been successfully used to de scrih je-some low-energy prop- 
erties of carbon nanotubes,E3E3o thoiLgli-^he situation 
in quantum wires remains controversialBaO 

The effect of boundary conditions on the low-energK 
properties of LL's was first considered several years ago.EZi 
Moreover, the interplay between boundary, finite-size, 
and thermal effects has beenpshewn to alter considerably 
the properties of the systemESo In particular, the zero- 
temperature critical behavior of the bulk always crosses 
over to a boundary dominated regime. These studies 
are important to explain the experimental results of tun- 



neling spectroscopy into one-dimensional systems. More 
recently, it has been proposed that one-dimensional sys- 
tems with gapless degrees of freedom and open boundary 
conditions form a new universality class of quantam crit- 
ical behavior called 'bounded Luttinger liquids'. E^l 

A particular kind of boundary effect emerges in the 
case of inhomogeneities. In general, an inhomogeneous 
LL is modeled by allowing the velocities of collective ex- 
citations Up and Ua and the correlation exponents Kp 
and Ka- to vary in space. The absence of conductance 
renormalization in long high- mobility GaAs wires, for in- 
stance, has been analyzed and explained in terms of an 
inhomogeneous LL model, where the Fermi Hquid leads 
are replacead-hy—a jion-interacting one-dimensional elec- 
tron gas.Ei[E3u3l2j'L£l Furthermore, LL's with different in- 
homogeneity profiles have also been used in the context of 
the fractional quantum Hall effect, to deS|G£ibe transitions 
between edge states at different fiilings,E3Ell or between 
an edge state and a Fermi liquidBa 

With an eye to practical applications as diodes or tran- 
sistors, researchers have recently ho j gun j to ifchlpipaite het- 
erojunctions of carbon nanotubesE3E3'E3llilc3'E2lEil which 
look especially promising. They happen to be another re- 
alization of an inhomogeneous one-dimensional system. 
Taking this idea one step further, we have been led to 
consider another kind of heterostructure: a superlattice. 
The effect of electronic correlations in superlattices was 
initiated through a one-dimensional Habhard-like model 
called a Hubbard superlattice (HSL) ,e3cjo consisting of 
a periodic arrangement where the Hubbard on-site repul- 
sion U is turned on and off in a repeated fashion. De- 
spite its simplicity, a number of remarkable features were 
found, in marked contrast with the otherwise homoge- 
neous system: local moment weight can be transferred 
from repulsive to free sites, spin density wave (SDW) 
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quasi-order is wiped out as a result of frustration, and 
strong SDW correlations (in a subset of sites) could set in 
above half-filling. Furthermore, the evolution of the local 
moment and of the charge gap, together with a strong- 
coupling analysis, showed that the electron density at 
which the system becomes a Mott insulator increases 
with the size of the free layer relative to the repulsive 
one. More recently, the possibility of a periodically mod- 
ulated hopping aJ-arbitrary filling and magnetization has 
been considered.c3 

In order to generahze the effects of a superlattice struc- 
ture in an interacting one-dimensional system, we con- 
sider here a general Luttinger liquid superlattice (LLSL), 
making at first no reference to the underlying microscopic 
details. We show how its low-energy properties bear 
strong resemblance to a conventional Luttinger liqui^ 
However, as in the case of bounded Luttinger liquids,E£l 
new effective parameters have to be introduced, which 
are the superlattice analogues of the spin and charge ve- 
locities and stiffnesses. These encode all the informa- 
tion necessary for a description of the low-energy sector. 
Moreover, these effective parameters turn out to mix the 
properties of the underlying sub-units in proportion to 
their spatial extent. This spatial averaging character- 
istic suggests the possibility of fine-tuning the physical 
properties by a careful selection of the superlattice mod- 
ulation, a feature which may prove useful in nano-device 
applications. We then consider specific realizations of the 
LLSL by analyzing in full detail a general HSL. We find 
a proliferation of phases, both metaUic and insulating. 
Surprisingly, the insulating phases often have no charge 
gap, because additional charge can be accommodated in 
the compressible sub-units. A partial account of these 
results has appeared in Ref. 

The paper is organized as follows: In Sec. II, we intro- 
duce the bosonic formulation of the Tomonaga-Luttinger 
model and our model. We obtain the effective charge 
and spin velocities, the correlation functions with the ef- 
fective exponents and the Drude weight for LL superlat- 
tices. The application of these results to various cases 
where the LL describes the low-energy sector of a Hub- 
bard model is analyzed in Sec. We close with the 
conclusions in Sec. IV. 



II. THE MODEL 



We briefiy review the general aspects of a homogeneous 
LL in order to set up the notation. The low-energy, large- 
distance behavior of a one-dimensional fermionic system 
with spin-indftpepdeijit-i^Bj^actions is described by the 
Hamiltonianll°lfiilJi^ 



H = H„ 



251 



(27ra)- 



■ j dxcos(V8$^), (1) 



where a is a short-distance cutoff, gi is the spin 
backward-scattering amplitude, and 



Hi, = dx 



2ttK, 



(2) 



with V = p and a for the charge and spin degrees of 
freedom, respectively. 
The phase fields are 



(3) 

and 

Here Pr{p) [fr(p)] are the Fourier components of the 
charge- (spin-) density operator for the right- (r — +) 
and left- (r — ~) branches of moving fermions. Introduc- 
ing the total number operators (measured with respect to 
the ground state) N^s for branch r and spin s, the total 
(charge and spin) number and current operators N^,Ji, 



,T)±(iV+^i+iV_,i)], (5) 



and 



= ;^[(^+^T - ^-,t) ± (N+.l - N-.i)], (6) 

where the upper and lower signs correspond to ly — p and 
a, respectively. 

The operators <i>,y and Hi, in Eqs. (|5) and (|^) 
obey Bose-like commutation relations: [^^(x),Ilf_i{y)] = 
iS,jf_i6{x — y). Consequently, at least for gi = 0, Eq. (0) 
describes independent long- wavelength oscillations of the 
charge and spin density, with Hnear dispersion relations 
Lu„{k) = Ui,\k\, {ui, is the velocity of elementary excita- 
tions) and the system is conducting. The only nontrivial 
interaction effects in (|]) come from the cosine term. How- 
ever, for repulsive SU(2) invariant interactions {gi > 0), 
this term is renormalized to zero in the long-wavelength 
limit, and at the fixed point one has K* — 1. The three 
remaining parameters in (|^) then completely determine 
the long-distance properties of the system; in particular, 
Kp determines the long-distance decay of all the correla- 
tion functions of the system. 

We now consider a LLSL, consisting of a repeated pat- 
tern of two different LL's with parameters ua,i^ K\^i, and 
sizes La (A = 1,2) perfectly connected (Fig. |l|). We use 
the adiabatic approximation, in which the scale of the 
inhomogeneity between the two liquids is much larger 
than the Fermi wavelength 2'jT/kp. Thus, the single- 
particle backscattering from the inhomogeneities can be 
neglected. Accordingly, the low-energy properties of this 



3 



1^2 U2v K which illustrate their duality under the replacement 

m ggggggggj V 2v ^ _^ (-^) ^ Substituting (§) into (|) yields 



.1 UivKiv 



Figure 1: Schematic representation of a Luttinger liquid su- 
perlattice. Here, ua,.^, Kx^i, and La are the velocities, inter 
action parameters and sizes of two Luttinger liquids (A=l,2) 



perlattice. Here, u\^^, K\^,y and L\ are the velocities, inter- dtt'^v ~ u^K^dx ( -^d^'^v ) — 0, (10) 



LLSL are described by generalizing the usual bosonized 
Hamiltonian of Eq. (m) as follows: 



and a similar equation for 0^. 



U—p.(7 



KJx) 



where the sum extends over separated charge- {v = p) 
and spin- (i/ = a) degrees of freedom, each of which 
with interaction- and layer-dependent parameters Ui,(x) 
and K^{x). For x on the first (second) 'layer' one has 
ify(x) = Ki^i,{K2,i,) and Ui,{x) = ui^„ ("2,1^)- 

The boson phase fields are related to the charge and 
spin densities, p and cr, through ^/2dx^u{x)/TT = while 
is such that dxQv is the momentum field conjugate to 
[$^(x),9ye^fy)] =i5^n5{x-y). Note that n^, (x) = 
dxQu{x) in Eqs. % and (§. 

The equations of motion for the fields and 0^ are 



dt'^i, = Uy{x)Ky{x)dxQu, 



K,{x) 



We now have to set up the matching equations at the 
interfaces between layers. The equations of motion lead 
to the continuity of and Qi, and their time deriva- 
tives. The right hand sides of Eqs. (§j|) yield, as addi- 

(7) tional conditions, the continuity of both {u„/K„) dx^v 
and UuK^dxQu at the contacts. Note that the conti- 
nuitjj— o£,|^ and 8,y guarantees that of the fermionic 
field.EM3€a Physically, these boundary conditions sim- 
ply encode the conservation of both charge and spin cur- 
rents ju = ^/2dt^u/iT (since we are neglecting Umklapp 
processes and backscattering of electrons with opposite 
spin). We stress that, under these conditions, these are 
the only universal requirements on the fields, irrespective 
of the actual interface potentials. 

The superlattice structure is incorporated into the so- 
lution of the equations of motion in a way completely 
analogous to the discussion of refiection and transmission 

(8) in the Kronig-Penney model. That is, we diagonalize the 
Hamiltonian (^ by expanding the phase fields in normal 

(9) modes 



^,{x,t) = _z^sign(p)^^[6_p,.e-»-'* 



Qu{x,t) = ^ y^ 



9p,y{x) 



(11) 
(12) 



where 6^ ^,, are boson creation operators (p > 0). The 
normal mode eigenfunctions 4>p.y{x) and eigenvalues LJp,i, 
satisfy 



i^l,ApAx)+'^i^K^dx{^dx(t)p,u]=Q, (13) LVp^^L2 .^^, u}p.uLi 
[obtained by taking ( pi] ) into (p^)], subject to the same 

boundary conditions at the contacts as before, with where A,y = 77,^ -I- 77"^ and 77^ = Kx^yjK^^v For p <t. 

<^p,y(x) replacing The eigenvalues are given by it /{Li-\-L2), the dispersion relation of the LLSL is linear. 



C0Sp[Li+L2) — cos[— )cos(— ^- ) 

U2,u 
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i.e., u}^{p) = \p\, with an effective velocity 

^1/ ~ j 5 



(15) 



where I = L2/L1; clearly, Ci, U2.y as I ^ 00, and 
Cv ui^^ as £ 0. Also, from Eq. (|lj) it follows that 
the spectrum of elementary excitations of a LLSL has 
bands and gaps, reflecting the superlattice structure. In 
this regard, it should be mentioned that, for a Luttinger 
liquid with a periodically modulated particle density, the 
presence of a plasmon gap was reported.E3 Here, we will 
focus only on the low energy properties of the LLSL. 

On the other hand, the zero mode functions (/)o,,y(x) 
and 9oy{x), satisfy 



K,{x) 

ixu = Ui,{x)K^{x)dx9QM{x), 



(16) 
(17) 



which follow from Eqs. (||) and While for the homo- 
geneous system one has 



In a LL, the ground state value of Tp measures the charge 
compressibility, whereas t„ is related to the spin suscep- 
tibility. Considering the LLSL zero modes [Eqs. (20) and 
(^l|) ] and the Hamiltonian (|^) we flnd that the superlat- 
tice compressibility is given by 



l + £ 



(25) 



where k\ ~ 2K\ p/T:u\^p is the compressibility of each 
layer. Clearly is nothing but an average of the indi- 
vidual compressibilities weighted by the layer lengths. 

Interactions in a one-dimensional system can enhance 
charge density or superconducting fluctuations depend- 
ing on whether they are repulsive or attractive. Let us 
then consider the correlation functions for the LLSL at 
T = 0. The asymptotic (i.e., for well separated x and y) 
behavior of the density-density correlation function is 



{n{x)n{y)) 



7r2 \x - yf 



^2!:(0(x)-</.(j,)) 



A. 



(26) 



and 



(18) 



(19) 



for the LLSL there will be, in general, an inhomogeneous 
periodic density proflle. As we will see, there is a ten- 
dency for the charge to accumulate more in the less in- 
teractive layer. Thus, the zero mode functions will reflect 
this inhomogeneityH Now, since each layer is a LL, (f)^^ 
and 6*0,1/ will vary in such a way that A0o,i/ — t^N\i, 
and A0o,i^ — ttJai^ across each layer A, with layer-speciflc 
number and current operators. We then obtain 



where 



[m - 1) 7rL2 



ttNx„x 
ttJxvX 



N2 



if A = 2, 



(20) 
(21) 

(22) 



with an analogous expression for Bm,x,i' obtained with 
the replacement of N\^„ by J\^v Here m — 1,2,3,.... 
labels the unit cell. Analogously, from Eqs. (^ and (||) 
we have 



7a,. = TTU^{x)Kt,{x)^, 



_ u^{x) Nx^ 
K,{x) La ■ 



(23) 
(24) 



where 



1 



K2 



f{Ki,,,K2 



(27) 



i + e 



) if X and y el, 
if {x,y) e (1,2), 



K2. 



^ "1,1/ 



K2.UU2U 



^ ) if X and y e 2, 



(28) 

and 4){x) = kpx—(j)o^p{x). The second and third terms on 
the right-hand side of Eq. (|2^ ) respectively correspond to 
the 2kF and Akp correlations in the homogeneous case. 
And, similarly to the homogeneous system, the former 
dominate over the latter for K* > i (see, however, Ref. 

The correlation functions for spin-spin, singlet (SS) 
and triplet (TS) superconducting pairing are given by 



(S(x).S(y)) 



Bi- 



B2- 



,2i{4,{x)-<P(y)) 



\x - y\ 



(29) 



Olsi^)Ossiy)) - (45o(^)Ot5o(2/)) (30) 



o^Ts^,{x)OTs^Ay) 



\x - y\ 

C2 

\x - y\ 



(31) 
(32) 



5 



where = l/i^a,,.) [Eq. reflecting the 

duahty properties (in the homogeneous Umit we have 
l/K^). One should note that the correlation func- 
tions depend not only on the difference x — y, but also 
on the actual positions x and y, through the zero mode 
functions. It is interesting to note that, even though 
we now have new effective coupling constants {K*,K^), 
the scaling laws between the exponents of the correlation 
functions are not broken by the superlattice structure. In 
other words, the replacement K* and K^, 

in the exponents of the correlation functions of the ho- 
mogeneous system yields the exponents given above for 
the superlattice. 

Finally, we discuss the conducting properties. Let us 
first consider a LLSL in the presence of a weak external 
space- and time-dependent electrostatic potential V{x, t), 
such that the electric field E{x,t) = —dxV{x,t). The 
interaction of the fermions with V{x,t) is described by a 
source term 



Hext ^ ~^ J dxp{x)V{x,t). 
Now the equation of motion for <i>p isBEliiii 



(33) 



Up (x) Kp (x) 



dx 



Up jx) 
Kp{x) 



$p(a;, t) = -eE{x, t). 



Defining the bosonic Green's function 

G{x,y,t) = -ie{t) {[^p{x,t),^p{yM) 
the nonlocal conductivity is given by 

cr{x,y,t) = -—dtG{x,y,t), 



(34) 
(35) 

(36) 



where go = e'^/h is the conductance quantum. First, we 
consider the usual order of limits, taking g — > before 
Lo ^ 0, which yields the Drude weight, appropriate for a 
situation of a uniform static electric field.E3 In this case 



aiq^O,Lo^O) = 2goCpK;S{u;), 



(37) 



which has the same form as for the homogeneous case,Eil 
but with the effective velocity and effective exponent re- 
placing the corresponding uniform quantities Up and Kp. 
Taking the limits in the reverse order yields the Landauer 
conductance, which corresponds to a situation where aa 
electric field is applied to a finite region of the sample.E^I 
In the LLSL we have 

aiq^0,u;^0)^2goK;S{q), (38) 

which is similar to the homogeneous case@ except that 
the effective exponent appears. Naturally, the conduc- 
tance renormalization of Eq. (QUA-usually hidden in the 
presence of Fermi Hquid leads EirE^O However, it should 
be accessible in AC measurementSj-if uj > Cp/L, the in- 
verse traversal time of the sample.E3 
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Figure 2: Schematic representation of Hubbard superlattices. 
In (a) the hopping is uniform but the interaction is (/a in the 
A sub-chain. In (b) the interaction is uniform, whereas the 
hopping can assume two values. 



III. HUBBARD SUPERLATTICES 

For the sake of illustrating the LLSL with a specific 
realization, we now--diS|ip^s-| a one-dimensional Hubbard 
superlattice (HSL) .a'C3l£dc3 We first consider a periodic 
arrangement of Li sites in which the on-site coupling is 
Ui > 0, followed by L2 others with on-site coupling U2 > 
0; the hopping parameter, t, is uniform, as shown in Fig. 
^(a). We subsequently consider the on-site interaction 
as being uniform but the hopping integrals as periodic: 
tl between Li sites, followed by t2 between L2 sites; see 
Fig. 1(b). 

Both cases above are contemplated if one writes the 
Hamiltonian as 

H = -'^U^i+i{cl^Ci+ia + H.c.) + ^^Uini^nii, (39) 



where, in standard notation, i runs over the sites of 
a one-dimensional lattice, cl^{cia) creates (annihilates) 
a fermion at site i in the spin state cr =t or | and 
n-ia = c|^Ci(j. It is important to notijee that the SL struc- 
ture breaks particle-hole symmetry.cll The homogeneous 
Hubbard model, in a grand-canonical ensemble descrip- 
tion, is invariant under a particle-hole transformation 



i-iy 



only when /i = U/2. In the superlat- 
tice case, a uniform chemical potential cannot ensure this 
symmetry throughout the whole system. Instead, under 
a particle-hole transformation the system is mapped onto 
a different one with a spatially modulated chemical po- 
tential. 

A weak coupling perturbation theory, similar to that 
for the homogeneous model can be used to show that 
Eq. (0) indeed describes the low energy and small mo- 
mentum sector of the discrete model of Eq. ( |39| ) in the 
limit of long layers; see the Appendix. Then, in Eq. (0) 
one has K^{x) = and u^{x) = ux^, for x on the 

layer A = 1,2, where and uai/ are the usual uni- 
form weak coupling LL parameters for each layer. It is 
by now well established that a LL description is appro- 
priate for the low-energy sector of the Hubbard model. 
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even in the strong coupling limit U c>o0 Now, each 
long Hubbard sub-chain is still a finite-sizedrXL, though 
connected to particle reservoirs at each end.Eil We there- 
fore make the quite reasonable assumption that the above 
LLSL description remains valid even in the strong cou- 
pling limit. With respect to magnetic properties, the 
superlattice structure (with repulsive interactions) does 
not break SU(2) symmetry, so that the inhomogeneous 
Ka- is still expected to renormalize to K„ 1. 

Because each sub-chain is an open LL, there will be 
a certain amount of charge redistribution between them, 
leading to a non-uniform charge profile. Let us first con- 
sider the special case of two layers only [with parameters 
{Ui,ti) and {U2,t2)] initially disconnected and with the 
same initial density n = N/L. In general, these two sub- 
systems will not have the same chemical potential. We 
then bring them in contact with each other, so that par- 
ticle exchange is allowed. Electrons will flow from one 
system to the other until their chemical potentials ex- 
actly match: 



fi{ti,Ui,ni) = /i(t2,C^2,"-2), 



(40) 



where /i and nx are the chemical potential and the equi- 
librium densities of each layer, respectively. This is just 
the condition for thermodynamic equilibrium. Naturally, 
conservation of total charge dictates that 



m +en2 ^n{l + i). 



(41) 



In order to determine ni and n2, we must solve simul- 
taneously Eqs. ( ^ and (41). The extension to the case 
of more than two layers leads to no modifications of the 
above equations and the charge profile will be periodic 
with the densities determined as above. 

The dependence of fi on the density n and on the in- 
teraction U can be obtained frouL-lhe exact solution of 
the homogeneous Hubbard model.EZl As a function of n, 
the chemical potential fi{t, U,n) increases monotonically 
and is discontinuous at half-filling, where it jumps from 
li-{t, U) to ^+(i, U) — U — f^-{t, U). Thus, the homoge- 
neous model is a Mott insulator at half-filling. ^_(t, j 
is the lower chemical potential at half-filling, given bj 



H-it,U) = 2t-U / 
Jo 



Ji{uj)duj 



(42) 



where Ji {lo) is a Bessel function. To increase the particle 
number above half-filling, we need to pay an energy given 

by 

= li+{t,U) - ^l-{t,U) = U - {t,U), (43) 

which is the quasiparticle gap. For later use, we also 
quote the chemical potential of the non-interacting case. 



fj,(t,0,n) 



-2tcos ( — j 



(44) 
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Figure 3: Phase diagram of a Hubbard superlattice with 
Ui = and ti — t2 = t for two layer length ratios: £ = 0.5 and 
2. For each £, there are two metallic phases and two insulating 
ones. The two dots locate {Uc,nc), where Uc/t — 3.2309 and 
nc^{2 + e)/{l+l). 



A. The Ui=0 case 

We first consider the case in which one of the layers 
is 'free' {Ui = 0) and take ti = t2 = t for simplicity. 
Figure ^ shows the phase diagram for £ — L2/L1 = 0.5 
and 2; the case £ — 1 has been discussed in Ref. For 
the sake of comparison, one should also keep in mind 
the phase diagram for the homogeneous LL, in which 
there is a single gapped (Mott) insulating phase for any 
non-zero repulsion at half-filling; upon either electron- 
or hole-doping the system becomes metahic. In what 
follows, we start with a qualitative discussion of the phase 
diagram, after which we provide the details of how the 
boundaries and special points are determined. 

In the case of a superlattice, while for n < 1 the sys- 
tem is always metallic, interesting metal-insulator tran- 
sitions have been found for n > 1, as displayed in Fig. 
^. Indeed, for a density n just above half-fihing, the sys- 
tem is still metallic, with more particles occupying the 
free layer than the repulsive one in order to decrease the 
overall electronic repulsion: One has ni > 1 and 712 < 1, 
as shown in Fig. As the density is increased for given 
£ and U2, electrons will be accommodated in both layers 
without affecting the metallic character; see Figs, js] and 
^. This will persist until the repulsive layer is half fihed 
(n2 = 1), when it becomes a Mott insulator. Recah that 
an insulating phase in one of the subsystems is signalled 
in Fig. Q by a horizontal plateau in the corresponding 
ni{fi) (i = 1, 2) plots. The system as a whole is therefore 
an insulator, since it can be thought of as a series ar- 
rangement of resistors. However, the unusual fact is the 
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Figure 4: Particle densities for the 1 = 2 Hubbard superlat- 
tice with U\ = 0, as functions of the chemical potential: n 
is the overall charge density (full lines), ni is the density at 
free sites (short-dashed curves), and 712 is the density at re- 
pulsive sites (long-dashed curves). Two cases are considered: 
(a) U2 = 2t< C/c, and (b) (/a = 4f > Uc. 



gapless nature of this insulating phase: charge can be 
accommodated in the free layer at no energy cost, since 
the system is compressible {dn/d^ 7^ 0) in this range of 
n; see Fig. |4[ 

As the density is further increased, the system re- 
sponds in two different ways, depending on whether U2 
is larger or smaller than Uc = 3.2309t (for all £); see Figs. 
I and |. If U2 < Uc [Fig. |(a)], the insulating state can 
only be sustained up to a limited amount of additional 
charge; that is, as long as it is energetically favorable to 
accommodate this extra charge in the free layer, while 
keeping n2 = 1. Further increase in n soon leads to an 
increase in the occupation of the repulsive layer (with 
2>ni>n2>l) and the system reenters an overall 
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Figure 5: The Mott-Hubbard gap at n = ric as a function of 
the on-site coupling for both the homogeneous model (dashed 
line) and the superlattice (full line). This behavior is the same 
for all £. 



metaUic phase. This metaUic character will be lost again 
for larger n, when the free layer becomes completely full 
{ni = 2, 1 < < 2), with the superlattice displaying 
insulating behavior. Again, this insulating phase is gap- 
less. 

If U2 > Uc [Fig. 1(b)], on the other hand, all added 
electrons will be accommodated in the free layer (1 < 
Til < 2,n2 = 1), so that the superlattice remains in the 
state of a gapless insulator. Further increase in the elec- 
tron density leads to the free layer becoming a band insu- 
lator (ni — 2) , while keeping the repulsive one pinned at 
half- fining; the density, Uc, at which this occurs depends 
on the aspect ratio, £, and is given by [c.f., Eq. ( pl| ) and 
Refs. |6|Q Uc = {2 + i) / {1 + i). Only at this special 
density does the superlattice become a Mott insulator, 
since it is incompressible {dn/dfi = 0); see Fig. 1(b). For 
n > Uc, the free layer remains completely full, so that all 
added electrons go to the repulsive layer; the superlattice 
behaves again as a gapless insulator. 

At this point it is worth commenting that the phase di- 
agram of Fig. ^ differs in two aspects from the one found 
for thin layers, obtained by means of Lanczos diagonal- 
izations: In Ref. 46 no gapless insulating phases were 
probed, and the insulating phase for n — Uc was found 
to extend down to any U2 > 0. The former difference is 
due to the fact that only gapped insulating phases were 
probed, while the latter can be traced back to finite-size 
effects. 

Let us now fill in the details on how the lines and spe- 
cial points of Fig. |^ are determined. The dotted horizon- 
tal line at n = Uc is obtained by setting both ni = 2 and 
n2 = 1. However, this condition can only be obtained if 
U2 > Uc, where Uc is determined implicitly by (see Fig. 
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Figure 6: The correlation exponent K* as a function of the 
density n for Ui = 0, [/2/t = 2 and t = 0.5, 2. 



M+(i,C/c) = M(i,0,ni=2) = 2i. 



(45) 



Since neither Eq. nor Eq. depends on £, this 
condition yields the same Uc ~ 3.23097t for any finite 
aspect ratio. 

Besides, for U2 < Uc, the system is always gapless. 
For U2 > Uc and n — Uc, the system shows a Mott- 
Hubbard gap given by the energy difference between the 
highest occupied state, which is the upper edge of the 
non-interacting band at 2t, and the lowest unoccupied 
level, which is the higher chemical potential of the half- 
filled Hubbard chain at fJ,+ {t, U2) 



Ai+(t, C/2) - 2i = f/2 - 2t - ^_ (i, C/2). 



(46) 



For the one-dimensional Hubbard model, one has A// ~ 
{%VtU /'K)exp{—2TTt/U) in weak coupling and /S.h oc U 
in strong coupling. For the HSL, we found that A5 is 
linear with U2 for large U2 and is always lower than the 
gap of the corresponding homogeneous system; see Fig.||. 

The two metallic phases are characterized by ni < 
2,712 < 1 (lower one) and ni < 2,ni > 1 (upper one). 
The metal-insulator transition (MIT) lines can therefore 
be obtained by means of Eqs. (pO|), (|4l| ), the Lieb-Wu 
chemical potential ^ (t, J72, '^^2) and (|44|). Therefore, in 
Fig. ^, (i) n' is the line in which the lower Hubbard band 
of the interacting sub-chain becomes fully occupied, (ii) 
n" is the one in which the upper Hubbard band starts to 
fill, and (iii) n'" is the line in which the non-interacting 
sub-chain fills up. Thus, 



M- {t, U2) 
{t, U2) 
2t 



,Uu2h^'^ 



(47) 
(48) 

(49) 



The LL description of Sec. H is only valid in the metal- 
lic regions of the phase diagram, where no gap is present 
in either the spin or the charge sectors. In these regions, 
we have 



1 + /\plVF/u2,p 



{tVFlu2,py 



1 -I- IvfIu2^p 



where vp is the Fermi velocity. When the insulating 
phase is approached from the lower metallic region (see 
Fig. H), Cp —> as a result of U2,p ^ in the interacting 
layer. In Fig. ||, we show the effective exponent if* as a 
function of the filling n. For any I, both metallic phases 
have 1/2 < K2^p < K* < 1 and the charge and spin cor- 
relation functions decay faster than in the homogeneous 
system. This is a direct consequence of the 'weighted 
average' character of the effective exponent K*. By the 
same token, for a given n on the lower metallic phase, K* 
decreases as £ increases. In the upper metallic phase, K* 
always tends to the non-interacting value of 1 as the up- 
per insulating region is approached; for the superlattice 
with larger £, K* reaches 1 at a lower overall density. 



B. The general case: U2 > Ui ^ 

We now consider a more general HSL, with dif- 
ferent non-vanishing coupling constants on each layer 
(Ui U2), while keeping the same hopping amplitude 
t throughout the lattice (Fig. §|(a)). Using once again 
the exact expression foi:^the chemical potential as a func- 
tion of both U and n,LZl we have determined the charge 
profile of the superlattice system. The charge tends to 
accumulate in the layer with the smaller coupling, which 
we choose to call layer 1. This is rather intuitive, since 
electrons decrease their mutual repulsion energy by fiow- 
ing into the less interacting layer. 

The phase diagram for this HSL is very rich. We ob- 
serve six different phases, three metallic (Mi, M2 and M3) 
and three insulating (Ii, I2 and I3), each characterized by 
its charge profile, as shown for three illustrative cases in 
Fig. ^ The topology of the phase diagram is thesame for 
any £ and the limiting cases Ui ^ (Section III A ) and 
Ui — U2 (homogeneous chain) are recovered. On each 
phase diagram of Fig. ^ there are five MIT lines, labeled 
by nj through ny , which are determined similarly to the 
case Ui = discussed before (see Table |). We get: 
line I: 



^i-{t,Ul) = ^l{t,U2, 



{l+£)n - 1 



line II: 



fl+{t,Ul) = fl t,U2, 



(l-l-^)n-l 



(50) 



(51) 
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Figure 7: Phase diagram for a ^ = 1 {ric = 1.5) Hubbard 
superlattice with Ui > 0, and three different values of U2. 
In each case there are three metallic (M) phases, and three 
gapless insulating (I) phases; gapped insulating phases appear 
as horizontal dashed lines. In (a), U2 = 3t and U* ~ 2.15t 
(see text); in (b), U2 = U* ^ 2.S9t, and Ua ^ O.Mt; and 
in (c), U2 = 16t, U* fa 3.02t, and Ua ~ 12. 2t. 



Table I: This table lists the various metallic (Mi , M2 and 
M3) and insulating phases (Ii, I2 and I3) of Fig. ^ with the 
corresponding sub-chain densities. The last column shows the 
nature of the transition lines between the phases (n/ through 
nv in Fig. |^. LHB A and UHB A respectively stand for lower 
Hubbard band and upper Hubbard band in layer A = 1, 2. 





Sub-chain densities 


Transition line 


Ml 


ni < 1, 722 < 1 








LHB 1 fills up (72/) 


Ii 


ni = 1, 722 < 1 








UHB 1 starts to fill (nj/) 


M2 


ni > 1, 712 < 1 








LHB 2 fills up {mil) 


I2 


721 > 1, 7l2 = 1 








UHB 2 starts to fill (72/^) 


Ms 


72i > 1, 712 > 1 








UHB 1 fills up (nv) 



line III: 



line IV: 



line V: 



fi{t,Ui,{i + e)n^e) = ^i^{t,U2); 



(52) 



fi{t,Ui,il + i)n~£) = ^^+{t,U2); (53) 



/2(i,[/i,ni = 2)===A._(<,C/2,^i±4^). (54) 



For Ui = 0, the lines njj j, njy , and ny deterpwne the 
phase diagram of Section [[IIA| [Eqs. (|^)-(g)]Jl3 

One of the consequences of a non-zero Ui is to push 
the lower metallic phase of Fig. ^ to smaller densities, as 
shown in Fig. (Mi). In addition to this phase, which 
spans all values of J7i < U2, there are two other metallic 
regions (M2 and M3). And in-between metallic phases, 
one finds insulating phases, one of which (Ii) is now sta- 
ble for n < 1, unlike the case for Ui = 0. These insu- 
lating phases have either n\ — 1, X — 1,2, or ni = 2 
(see Table |). Once again, there is a 'division of labor' 
between the two types of sub-chains: while one is gapped 
(Mott) or completely filled (band), being responsible for 
the insulating behavior of the system, the other remains 
gapless and so does the system as a whole. 

Figure 0(a) shows the phase diagram for U2 — 3t < Uc 
(Uc is the same as for the case Ui = 0). The HSL has a 
gap at the density n = 1 iov Ui > U* 2.145608i; this 
n = 1 line separates the Ii (i.e., ni = 1, 722 < 1) and 
the I2 (ni > 1, 722 = 1) gapless insulating phases. For 
Ui < U* one goes through a sequence of MIT's, in which 
all insulating phases are gapless. 

In Fig. 0(b), we show the phase diagram for U2 — 
4t > Uc- As the overall density is increased from 1 in 
the interval Ua < Ui < U* , where Ua « 0.6433i and 
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Figure 8: Parameters U* and Ua as functions of U2. The two 
curves meet at f/2 = 6. 2526 It. 



Figure 9: The Mott-Hubbard gaps AJ and As. a at the den- 
sities n — 1 and ric = {2 + £)/{! + £), respectively. Here, 
U2 = 16t and 1=1. 



U* « 2.39149t, the system goes through a sequence of 
MIT's without ever being gapped. However, for Ui < Ua, 
the intermediate I2 {ni > 1, n2 — 1) and the I3 (rti = 2, 
712 > 1) gapless insulating phases are separated by the 
dashed hne at the density Uc = {2 + £) /{! + £), where the 
system is fully gapped. Similarly, another gap appears 
at the density n = 1 for Ui > U* , which again separates 
gapless insulating phases Ii (ni = 1, n2 < 1) and I2 

(ni > 1, 77,2 = !)■ 

For U2 = 16t > Uc [Fig. 0(c)] and U* < Ui < Ua (now 
U* = 3.01509t and Ua = 12.1724i) the system is metallic 
only below line /, which approaches n = 1 for large Ui; 
also, gapped behavior is again observed at densities n — 1 
and Uc, with all other insulating phases being gapless. 
For each of the regions Ui < U* and Ui > Ua, & 'tipped' 
metallic phase is observed. 

The above discussion indicates that there are special 
values of Ui, U* and Ua, which respectively represent 
the 'tip' positions of the low- and high-density metallic 
phases. Their dependence on U2 can be extracted from 
the solutions of 

fi+{t,U2)^fi{t,Ua,ni^2)^2t + Ua (55) 

and of 

/i+(i,f/*) = /i+(t,/72), (56) 

and are shown in Fig. ^. It should be noted that these 
values are independent of the aspect ratio £. As Fig. ^ re- 
veals, one should not be misled by the different horizontal 
scales in Fig. |^: the low-density tip does not recede as 
U2 increases, since U* actually increases monotonically 
with U2, saturating at Uc as U2 ^ 00. On the other 
hand, Fig. ^ shows that Ua is only defined above a cer- 
tain threshold, U2 = Uc, reflecting the fact that when the 



coupling in layer 2 is small, the situation ni = 2, 712 = 1 
is never realized; above Uc, Ua increases linearly with U2. 

According to our previous analyses, these two curves 
(which intersect at U2 = 6.25261i) deflne regions in the 
(Ui ,U2) plane characterized by the number of gaps in the 
sub-units for appropriate fillings, as specified in Fig. ^. 

Similarly to the case Ui — Q, the gaps at the densities 
n = 1 and n = Uc are given, respectively, by 

= ^I+{t,Ul)-^l-{t,U2), (57) 

As,a = fi+{t,U2)~2t-Ui, (58) 

and, again, they do not depend on £. The gaps and 
As.a, for U2 = IQt, are shown in Fig. || as functions of Ui. 
The gap at n = 1 [n = Uc] increases [decreases] linearly 
with Ui and vanishes for Ui < U* [Ui > Ua]. 

For the Hubbard model w^h repulsive interactions we 
have Ua < vf and Up > vpl^ For U2 — 4,t and Ui = 2t, 
the effective charge and spin velocities for the £ = 1 one- 
dimensional Hubbard superlattice are shown in Fig. |l^ 
as functions of n. The effective charge velocity (full line 
in Fig. ^ vanishes upon approaching the insulating 
regions as a result of the vanishing charge velocities of 
the individual sub-chains ux,p 0. Thus, Cp shows a 
re-entrant behavior as a function of n (cf. Fig. |^). As 
in the homogeneous case, the effective spin velocity is 
always smaller than the Fermi velocity and only vanishes 
in the upper insulating phase (dashed line in Fig. ^o|). 
The different behaviors of Cp and Ca can be traced back 
to the the fact that if* is sensitive to the superlattice 
structure, while if* — 1, since Ki„ — K2a = 1 as a 
result of the SU(2) symmetry being preserved. 

The preservation of SU(2) symmetry also leads to 
K„ = K* = I. Thus, from Eqs. (|§ and (||), the 
density-density and spin-spin correlation functions for 
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Figure 10: The effective charge and spin velocities Cp (full 
line) and Ca (dashed line) for a Hubbard superlattice as a 
function of n for Ui = 2t, U2 = 4t and £ — 1. 
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Figure 11: The correlation exponent Kp as a function of the 
band filling n for ^ = 1 and several values of the coupling 
constants. 



the HSL are dominated by {O^O) ~ |a; — yp^"^'" . These 
terms correspond to 2fcir-CDW and 2fci;'-SDW in the ho- 
mogeneous system. Here, i^2p < K* < Kip and the 
density-density and spin-spin correlation functions for 
the HSL decay faster (slower) than for a homogeneous 
system with U = U2 {U ^ Ui). Similarly, pairing cor- 
relation functions are {O^O) ^ \x — y\ ^ In spite 
of the presence of effective exponents K* and Kp, the 
condition for superconducting quasi-long range order is 
again K* > 1, analogous to the homogeneous case; this 
condition, nonetheless, remains unsatisfied. 



In Fig. |Tl|, the correlation exponent K* of the HSL is 
shown as a function of band fiUing, for different £ = 1 
superlattices: HSL-1 with U2 = 4t and Ui = 2t; HSL-2 
with U2 = 16t and Ui = 2t; HSL-3 with U2 = 16i and 
Ui = 8t. For any £, all metallic phases are characterized 
by 1/2 < < 1. We note that HSL-1 has three metal- 
lic phases, HSL-2 has two metallic phases and HSL-3 has 
only one metallic phase. On the low density sid e, K* ap- 
proaches 1/2 in contrast to the case Ui = (Sec. Ill A ), in 
which K* remains between 1/2 and 1. From Eg. (^7]) one 
sees that K* interpolates monotonically between Kip and 
K2p as £ is varied from to 00, highlighting the possibil- 
ity of a continuous 'modulation' of a physical parameter 
through the tuning of the superlattice structure. 



C. Two different hoppings: t2 > ti and U\ = U > 

We now consider two Hubbard chains arranged period- 
ically with the same couplingJ/A = U > 0, but different 
hoppings t2 > h (Fig. 1(b) ).EI 

Initially, the charge tends to accumulate in the layer 
with larger hopping (layer 2), because its chemical poten- 
tial is the smallest. Eventually, their chemical potentials 
become equal at the special density n*, determined by 
n{tx,U,n*) = 0. Then, for n > n* , the charge flow is 
reversed and proceeds from layer 2 to layer 1. n* is in- 
dependent of r = t2/ti and £, and decreases with U (see 
Fig. HI). 

It is interesting to plot a phase diagram in terms of 
the density and the ratio between the two hopping ampli- 
tudes, r = ^2/^1. We then identify seven different phases, 
three metallic (Mi, M2 and M3) and four insulating (Ii, 
I2, 13 and I4), as shown in Fig. ^ for [/ = Ati and U = &ti 
with £ = 1 and listed in Table I^ The value of n* lies 
within the Mi phase. We mention that several differ- 
ent insulating and metallic phases were also [found in a 
p-merized Hubbard chain in a magnetic fleld.ca 

Following the same reasonings as before, the lines in 
the phase diagram in Fig. |2| are given by: 
line I: 



fj.(ti,U,ni = 0) = ^ t2,U, 



n'j{l+£) 



line II: 



^i-itl,u)^^[ t2,u.. 



n'jj{l+£)~l 



line III: 



fl+iti,U) = fl t2,U, 



n'jjjil + £)-l 



line IV: 



line V: 



fi{ti,U,n'jy{l +£)-£) ^ pL-{t2,U), 



fi{ti,U,n'y{l +£)-£) ^ fi+{t2,U), 



(59) 
(60) 

(61) 

(62) 
(63) 
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Figure 12: Phase diagram for a Hubbard superlattice with 
different hoppings t\ and t2, in terms of filling n and hopping 
ratio tilt\\ the on-site coupling is homogeneous and set to 
U = Ati in (a) and U = 8ti in (b). There are three metallic 
(M) and four insulating (I) phases {£ — 1 and ric = {2+i)/{l + 
i) ). In (a), r* = 1.94847, = 2.26984 and rt = 4.4227 and in 
(b), r* = 3.818576, = 1.57735 and n = 6.02322. Besides, 
the values of n* are (see text for definition) (a) n* = 0.670511 
and (b) n* = 0.594067. 



line VI: 



/z(il,C/,ni=2)=/i t2,f/: 



n'yjil + £) 



(64) 



12 is 



Again, the topology of the phase diagrams in Fig. 
the same for any £. 

At small densities (Ii phase), charge accumulates in 
layer 2 while layer 1 is empty (ni = 0); the system is 
therefore a gapless insulator. 

As the density increases, layer 1 only starts being filled 
at n = n'j{r), determined by Eq. (|59|), which locates a 
transition to a metallic state (Mi); see Fig. |lj. Further 
increase in the overall density leads to an increase in both 
rii and n2. When layer 1 becomes half-filled, which oc- 
curs at n = n'jj{r) as determined from Eq. (pQ), the 



Table II: This table lists the various metallic fMi . M2 and 
M3) and insulating phases (Ii, I2, I3 and I4) of Fig. |l^with the 
corresponding sub-chain densities. The last column shows the 
nature of the transition lines between the phases (n'j through 
n'vi in Fig. |l2|). LHB A and UHB A respectively stand for 
lower Hubbard band and upper Hubbard band in layer A — 
1,2. 





Sub-chain densities 


Transition line 


Il 


7^1 = 0, n2 < 1 




V 




LHB 1 starts to fill (n'j) 


Ml 


ni < 1, n2 < 1 








LHB 1 fills up {n'ji) 


I2 


ni = 1, n2 < 1 








UHB 1 starts to fill (n'uj) 


M2 


ni > 1, n2 < 1 








LHB 2 fills up {n'lv) 


I3 


ni > 1, n2 = 1 








UHB 2 starts to fill (71^) 


M3 


ni > 1, 712 > 1 








UHB 1 fills up {n'vj) 


I4 


ni = 2, 712 > 1 





system reenters a gapless insulating state (I2). If r < r*, 
where 



^i-{h,u)^ 



(65) 



upon increasing the density the system goes through a 
gapped phase &t n = 1. The dependences of r* with U , 
and of the gap at ti = 1 , 



Ai+(ti,t/)-Ai-(t2,C/), 



(66) 



with r, are shown in Figs. ^ and |l^, respectively; note 
that r*{U = Atx) = 1.94847 and r*{U = Sh) = 3.818576. 
By contrast, if r > r* , the system enters a metallic phase 
(M2) bounded by njjj{r), and n'jy, given by Eqs. ( |6l] ) 
and (H). 

When increasing the density above half fiUing, the 
sequence of phases depends crucially on whether r* is 
smaller or larger than 



M+(i2,C/) 



(67) 



which, according to Fig. occurs when U < U ^ 
4.4191^1 or when U > U, respectively. 

Let us first consider U < U, which is the situation 
of Fig. ^(a). If 7- < r* < Tc, one goes through two 
transitions as n increases: I3 M3 at n'y [see Eq. (|63|)], 
and M3 I4 at n'yj [Eq. (|6|)]. If r* < r < r^, the 
sequence is M2-I3-M3-I4, until the lattice is completely 
filled. Another regime is determined by 



n{tuU,ni=2) 
M-(<2,C/) - 



(68) 
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Figure 13: Parameters r*, Tc and Ti as functions of (7/ti. 
Here, 17 = 4.4191ti. 




Figure 14: The Mott-Hubbard gaps A* and Ac.r at densities 
n = 1 and ric — (2 + £)/(!+ £), respectively. Here, U — 8ti 
and £ = 1. 



whose dependence on U jtx is also shown in Fig. |l^. If 
Tc < r < n, the system goes from a metallic (M2) to 
a gapless insulating phase (I3), and then, at n = ric, 
another Mott-Hubbard gap opens, which is given by 



Ai+(t2,C/)-M(ii,C^,«i = 2). 



(69) 



For a fixed ratio C//ii, A^r behaves as shown in Fig. 
[1^ . Above ric, the gapless insulating state I4 is again 
stabilized. It should also be noted that both gaps (A* 
and Ac,r) display universal behavior in the sense that 
they do not depend on I. Note also that Eqs. (ca), (67), 



0.5 





— r=2, U/t^=4 




- r=2, U/t^=8 




-- r=4, U/t^=4 


X 














1 , 


'1, 



0.0 0.5 



1.0 
n 



1.5 2.0 



Figure 15: The effective correlation exponent K'p {£ = 1) as a 
function of the band filling n for different values of U/ti and 
r. 



and ( pq ) do not depend on £, so that U is also universal. 

We now consider U > U, an example of which is shown 
in Fig. |2|(b). For r < rc < r* , one finds the same se- 
quence I3-M3-I4, with all insulating phases being gapless. 
If Tc < r < r*, a gapped insulating phase is crossed at 
n = ric- Similarly, for r > r* one goes from a metallic to 
a gapless insulating phase (M2 ^ I3), and again crossing 
the Mott-Hubbard phase at Uc- 

The effective charge and spin velocities are given by 



^ I J' 



(70) 



which vanish for n < n* and are smaller than the veloci- 
ties of the homogeneous system {d, < u^). Furthermore, 
Cp displays re-entrant behavior as a function of n. 
Finally, the effective interaction parameter if* is 



K* = 
p 



1 -I- A„£rui^i,/u2,u + (f^ui,i//u2,i/) 



1 

Ki,, 



7 1 

■K2,, 



(71) 



In Fig. |l5|, K* is shown as a function of band fiUing, 
for different coupHngs in superlattices with £ = 1: HSL- 
A with U = 4ti and r = 2; HSL-B with U = 8ti and 
r = 2; HSL-C with U = Ati and r = 4. Note that 
1/2 < if* < 1 for any £ in the metaUic phases. The 
various cases depicted in Fig. |l^ show three (A), one (B) 
and two (C) metallic phases. In the homogeneous Hub- 
bard chain, the density-density and spin-spin correlation 
functions decay faster when the hopping increases, since 
Kp increases with the ratio t/U. The effective correlation 
exponent of HSL-C is larger than in HSL-A (see Fig. pIs] ), 
because of the larger hopping amplitude of sub-chain 2 
in HSL-C and the 'averaging' nature of K*. 

We should stress that in the homogeneous system, the 
Luttinger Liquid description breaks down at half-filling. 
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when a gap opens in the charge (though not in the spin) 
sector. In the superlattice, this breakdown occurs in the 
insulating phases, as a result either of Umklapp processes 
(Mott gap, lower phase of Fig. ^, phases Ii and I2 of 
Fig. and phases I2 and I3 of Fig. ^2|), or of a band in 
one of the sub-lattices becoming completely full or empty 
(upper phase of Fig. ^, phase I3 of Fig. ^ and phases Ii 
and I4 of Fig.[l|), 



IV. CONCLUSIONS 

We have discussed in full generality the properties of 
Luttinger Hquid superlattices. We have seen how most 
features of a conventional Luttinger liquid description 
survive in the superlattice structure. In particular, a few 
effective parameters, the spin and charge velocities (cp 
and Cc) and the stiffnesses {K* and K^) are all that is 
required for a complete description of the low-energy sec- 
tor. These turn out to be combinations of the LL param- 
eters of the superlattice sub-units combined in proportion 
to their spatial extent. As we have stressed in the Intro- 
duction, this opens the way for possible 'engineering' of 
Luttinger liquids. 

This framework was applied to the study of the gen- 
eral phase diagram of Hubbard superlattices. It was then 
illustrated how one can tune between different phases by 
an appropriate choice of superlattice modulation. It was 
found that the superlattice displays a variety of metallic 
and insulating phases, the most prominent feature being 
the appearance of gapless insulating phases, as a result 
of the one-dimensional character of the system; gapped 
insulating phases were also found at some special densi- 
ties. 

Single- wall metallic carbon nanotubes (SWMN's) seem 
to provide a promising route towards realizing these 
LLSL's. Indeed, notwithstanding the fact that SWMN's 
are, in general, described by aJftss. simplistic model (pos- 
sibly even with more branchestaEa) , the LL coupling con- 
stant depends on its (true) aspect ratio throughEJ 

K^=ll + J^ln—\ (72) 

where k is the dielectric constant, and L and R are, 
respectively, the nanotube length and radius; typically 
one has Kp ~ 0.2-0.3. More recently, the growth 
of intramolecular junctions of SWMN's with different 
radii has been achieved with the introduction of a 
pentagpj|--amj_aJtie»tagon into the hexagonal carbon 
lattice,BElJ!3Jii'BEy so that the fabrication of a su- 
perlattice made up of SWMN's with different coupling 
constants has become a concrete possibility. 

We therefore expect the phase diagram of this 'nan- 
otube array' to share several features with the general 
Hubbard superlattice. This is because the only ingre- 
dients that enter into the phase determination are the 



thermodynamic equilibrium condition and charge con- 
servation. In the case of a Luttinger Hquid these can be 
easily written down if one knows how the LL parameters 
depend on the density 

Thus, the sequence of insulating and metallic phases that 
we have found in Hubbard superlattices should be present 
in other systems as well, as they will reflect the phase 
diagram of the sub-units. We hope this rich variety of 
behaviors will stimulate further experimental work along 
the Hues of carefully controlled nanotube arrays. 
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Appendix A: WEAK COUPLING BOSONIZATION 
OF A HUBBARD SUPERLATTICE 

Here we consider a Hubbard superlattice in weak cou- 
pling and show that it is possible to describe the low- 
energy properties in terms of a Luttinger liquid superlat- 
tice. The Hamiltonian of a Hubbard superlattice is 

H = -iI](*L*j+i.- + /i-C.) + IZ^:'"j.T^ja 

We focus on the low energy modeSpUear the Fermi surface, 
so that each fermion is written aalj 

^-j- , « e-''=^^'^*_ + e*'=^J'''*+j,^, (A2) 

where a is the lattice parameter, and the subscripts -I- and 
— respectively denote right and left movers. The kinetic 
energy part is then linearized as in the homogeneous case 

Ho = -t^(*,V*.+i,. + ft.c.)-ME*L*^.- 

(A3) 

The fermiooLc fields are given in terms of the bosonic 
ones, <&±,CT, aalj 

f±.<,(x) = -^C/i,,eT2.0F*±.„^ 
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Here-ct-is a cutoff parameter and U±,a- is the Klein 
factor .EjM Thus we get 



where = (<I>| ± ^i)/V^ and 



(A4) 



(A5) 



We now work out the low energy par t of the on-site 
Hubbard interaction. Again, we use Eq.(A2) to get 



i 



(A6) 



where : ... : denotes normal ordering,E^ J±,j,cr ='■ 
J- jr^'i j> -^dx'^±,(y, and the Umklapp terms have 
been neglected. Then 

H,nt ^ a [ dx C/(a;)-(a,$T)(5,$x) 



+ a dx U{x) 



(A7) 



where h.c. stands for hermitian conjugate. In terms of 
charge (p) and spin (cr) fields we have 



dx , 



a / dxU{x) 



2(7ra)2 



cos(V87r<I>^), (A8) 



the last term corresponding to the spin backscattering 
interaction, which is irrelevant in the RG sense. Finally, 
the low energy Hamiltonian for the Hubbard superlattice 
is 



H = 



1+ 

1 - 



U{x) 



ttvf 
U(xl 

TTVF 



(A9) 



This has the same form as Eq. (^, which describes the 
Luttinger liquid superlattice. 
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